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The interaction between light and matter is 
fundamental to developments in quantum optics 
and information. Over recent years enormous 
progress has been made in controlling the in- 
terface between light and single emitters includ- 
ing ions [1], atoms [2 ) 3], molecules [4], quan- 
tum dots [5], and ensembles |6j. For many sys- 
tems, inter-particle interactions are typically neg- 
ligible. However, if the emitters are separated 
by less than the emission wavelength, A, resonant 
dipole— dipole interactions modify the radiative 
decay rate |7J and induce a splitting or shift of the 
resonance [8j . Here we map out the transition be- 
tween individual dipoles and a strongly interact- 
ing ensemble by increasing the density of atoms 
confined in a layer with thickness much less than 
A. We find two surprising results: whereas for 
a non-interacting ensemble the opacity increases 
linearly with atomic density, for an interacting 
ensemble the opacity saturates, i.e., a thin dipo- 
lar layer never becomes opaque regardless of how 
many scatterers are added. Secondly, the relative 
phase of the dipoles produces an abrupt change 
in the optical transmission around the thickness 
A/4. 

The requirement for strong dipole-dipole interactions 
— that the atoms must separated by less than A — places 
stringent demands on the spatial extent of the ensemble. 
For atomic gases considerable progress has been made 
in the ability to confine atoms in thin layers [9] [10] , op- 
tical fibers pTJ [12] and waveguides [13] . Alternatively, 
it has been recognised that the spatial constraints can 
be relaxed by exploiting the long wavelength dipoles as- 
sociated with Rydberg transitions p~4j [15] where strong 
dipole-dipole interactions between Rydberg atoms give 
rise to a large optical non-linearity p~6j [T7]. However, 
micron-scale spatial confinement is still required [18] . As 
in the two-level case, in Rydberg systems dephasing can 
also play a role in the emission process [19] [20] . 

The underlying mechanism of light scattering is the 
interference between the incident field and the field 
produced by induced oscillatory dipoles. For indepen- 
dent dipoles, on resonance, the destructive interference 
in the forward direction is equivalent to removing a 
cross-sectional area of the incident beam equivalent to 



<Jo = 67rA 2 per atom, where A = X/(2ir) is the reduced 
wavelength of the dipolar radiation, corresponding to the 
effective length scale of the light-matter interaction. If 
the angular frequency uj of the incident field is detuned 
from the resonance frequency c^o by A = u — c^o, the 
optical cross section becomes 

a i + 4(A/r) 2 ' 1 ; 

where T is the full width at half-maximum of the res- 
onance. The frequency dependence of the opacity of a 
single dipole is illustrated in Fig. lb. For a dilute dipolar 
layer of thickness Sz and TV dipoles per unit volume, the 
opacity (optical depth) is given by 

OD = NaSz . (2) 

The optical depth increases linearly with iV, as shown in 
Fig. lb. However, when the density reaches a value where 
the dipole spacing TV -1 / 3 ~ A, the resonance is broad- 
ened and the efficiency of each dipole is reduced due to 
dipole-dipole interactions (Fig. lc-d). This effect is anal- 
ogous to dipole blockade in Rydberg gases [14] but for a 
thermal gas of two-level atoms it appears as an effective 
increase in the linewidth r = To + Tdd, where Tdd = @N 
is typically referred to as self-broadening [2TJ [22] and f3 is 
the self-broadening coefficient. Self-broadening results in 
a decrease in the resonant cross-section ao ~ 1/ (NX) (see 
Supplementary Information). In the high density limit 
the increased scattering due to adding more dipoles is 
exactly cancelled by the higher damping rate due to res- 
onant dipole-dipole interactions, leading to the surpris- 
ing result of a constant broadband opacity (see Fig. lc-d) 
given by (see Supplementary Information), 

OD ~ ^ . (3) 
A 

As an interesting comparison, if we take a monolayer 
of atoms (the thickness Sz is of the order of the Bohr 
radius, ao) and take the transition wavelength as A ~ 
ao/a, where a is the fine structure constant, we obtain 
an opacity OD ~ a (see Supplementary Information), 
very similar to the result of na recently demonstrated 
for graphene [23] . 

To confirm the saturation of opacity we measure the 
optical transmission of a single dipolar layer, where the 
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Figure 1. Evolution of opacity with density. Dipoles are represented in (a) by semitransparent discs, whose radii are the 
dipole-dipole interaction range A. At low density the discs do not overlap on average, and hence the interactions are negligible, 
whereas at high density the interactions dominate, (b) With no dipole-dipole interactions the lineshape is a Lorentzian with 
constant width and the sample quickly becomes optically thick (black), (c) When dipole-dipole interactions are considered, 
the line broadens linearly with density, eventually leading to a constant broadband value at high density, (d) Experimental 
data for cell thickness 5z — 390 nm showing transmission over the Rb D2 line, where there are now four separate absorption 
lines owing to the hyperfine splitting of the ground states of 85 Rb and 87 Rb. Opacity is represented by the colour map, which 
goes from completely transparent (white) through to completely opaque (black). 



Photodiodef^x 



Density 
Filter 




30 nm 



Experimental Cell 

^ Lens ^ , _,. 

^ 2 |jm^* 




Figure 2. Experimental setup and cell, (a) Schematic of 
optical setup. A laser is scanned across the D2 atomic reso- 
nance line in Rb (wavelength 780 nm). The light is filtered 
to low power to avoid optical pumping effects and focussed to 
a 30 /im spot size inside the cell, leading to a local thickness 
variation that is limited by the surface flatness of the win- 
dows. The transmission is then recorded on a photodiode. 
(c) The 'Newton's Rings' interference pattern is observed as 
the gap between the two sapphire windows changes thickness 
(b), due to the curvature of one of the windows, with a ra- 
dius R > 100 m. At the centre of these rings the cell has a 
thickness Sz ~ 30 nm. 



layer thickness must be less than or of the order of A. In 
our experiment we confine a gas of Rb vapour within a 
nanocell formed by two super-polished sapphire surfaces 
with a separation between 30 nm and 2 /im, similar to 
those in refs. [9j [10] • To achieve such small separations 
the plates have a surface flatness of less than 3 nm and a 



radius of curvature R > 100 m. After several iterations of 
cell design it has been possible to fabricate a nanocell of 
exceptional quality that allows the scaling of atom-light 
interactions in a dipolar layer to be elucidated. 

The experimental set up and the nanocell are shown in 
Fig. 2. A narrow band laser with a linewidth of less than 
1 MHz is frequency tuned in the vicinity of the D2 res- 
onance lines in Rb (A = 780 nm), and the transmission 
through the cell is recorded at different temperatures. By 
changing the cell temperature between 20° C and 350° C 
we can vary the number density between NX 3 <C 1 and 
A^A 3 > 100 . In doing so we tune the energy of dipole- 
dipole interactions (Vdd N) over 6 orders of magni- 
tude. The measured transmission as a function of laser 
detuning is shown in Fig. Id. The similarity with the 
interacting case in the high density limit (Fig. lc) is im- 
mediately apparent. At low temperature the spectrum 
consists of four main Dicke narrowed resonance lines cor- 
responding to the two isotopes, 85 Rb and 87 Rb, that are 
both split by the ground state hyperfine interaction (the 
excited state hyperfine splitting is also resolved in the 
case of 87 Rb). As the number density is increased this 
structure is completely lost and one observes a broadband 
absorption that saturates at a finite non-zero value. 

To make a more quantative comparision, in Fig. 
3 we plot transmission spectra for a cell length of 
A/2 = 390 nm together with the predictions of a the- 
oretical model that includes the effect of dipole-dipole 
interactions. The theoretical calculation is based on a 
weak probe absorption model [24 including the effects 
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Figure 3. Transmission spectra - experiment and theory. Ex- 
perimental (a) and theoretical (b) transmission spectra at 
thickness Sz — 390 nm for measured resevoir temperatures of 
130 °C (black), 160 °C (purple), 220 °C (blue), 250 °C (green), 
280 °C (yellow), 310 °C (orange), and 330 °C (red). Theory 
spectra are based on fitting the experimental data, using a 
Marquardt-Levenberg method (see [30]). Zero on the detun- 
ing axis represents the weighted line centre of the D2 line. 

of Dicke narrowing [25j [26] , cavity effects between the 
cell walls [27 and dipole-dipole interactions (see [28] and 
references therein). Although this is thicker than a sin- 
gle layer, the saturation of opacity can still be observed. 
Additional broadening and shifts due to van der Waals 
(vdW) atom-surface interactions are also present for 
Sz < A [29]. We restrict our investigation to Sz > 90 nm 
to avoid vdW atom-surface interactions strongly influ- 
encing the resonance lines, such that in the high den- 
sity limit of interest to this work dipole-dipole broad- 
ening dominates. For a random distribution of two-level 
dipoles the shift of the resonance (Lorentz shift) is smaller 
than the broadening [22] and is neglected. The effect of 
the Lorentz shift is apparent in Fig. 3a, but does not 
change the on resonance opacity significantly. Overall, 
Fig. 3 suggests that the model provides an excellent de- 
scription of the measured transmission. 

To illustrate the saturation effect, we plot in Fig. [4] 
the scaled opacity NaX = Im[x], where x 1S the sus- 
ceptibility of the medium (see Supplementary Informa- 
tion), as a function of number density. We take the 
opacity at the unshifted resonance position of the 85 Rb 
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Figure 4. Saturation of opacity due to dipole-dipole interac- 
tions. Without dipole-dipole interactions (dashed line) the 
opacity increases linearly with density. When interactions 
are included into the standard theory (upper solid line) the 
scaled opacity for the D 2 line is predicted to saturate at a 
value NaX = y/2. However for thicknesses Sz < A/4, we ob- 
serve a change in the behaviour which doubles the broadening 
coefficient (lower solid line). A constant additional broad- 
ening of r a d = 27T x 200 MHz has been added to account 
for atom-surface effects, which does not affect the value at 
which the opacity saturates at high density. Experimental 
data: Sz = 90 nm (magneta pentagons); Sz = 110 nm (red 
squares); Sz = 140 nm (green circles); Sz = 180 nm (blue tri- 
angles); Sz = 220 nm (yellow diamonds); Sz = 250 nm (cyan 
hexagons); Sz = 390 nm (purple stars). 

F = 3 — » F r = 4 hyperfine transition, though in principle 
one can choose any detuning and achieve similar results. 
For layer thicknesses A/4 < Sz < A, the data match our 
theoretical model including interactions extremely well. 
Also shown is the theoretical prediction without interac- 
tions, where the opacity increases linearly with density 
and does not saturate. For the D2 line, the saturat- 
ing value is NaX = \/2 following standard theory (see 
Supplementary Information). However for thicknesses 
Sz < A/4 we observe a deviation from this theory con- 
sistent with a dipole-dipole interaction strength j3 that 
is twice the expected value [28]. This arises because for 
Sz < A/4 atoms moving in opposite directions no longer 
dephase along the light propagation axis (see Supplemen- 
tary Information). Remarkably there is a sharp transi- 
tion between the two regimes at Sz = A/4 = 195 nm, as 
shown in Fig. [5] where we plot the effective strength of 
the averaged dipole-dipole interaction, scaled in units of 
the self-broadening parameter f3. 

In summary, we have demonstrated that adding more 
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Figure 5. Dipole-dipole interaction strength. Interaction 
strength f3 is plotted as a function of layer thickness. Values 
and errors are calculated from fitting data of the form shown 
in Fig. fusing a Marquardt-Levenberg method. We calculate 
the average interaction strength for Sz < A/4 (red circles) to 
be (2.2 ± 0.2)/3, whereas for 5z > A/4 (blue squares), we find 
(0.9±0.1)/3. 



scatterers in a thin layer does not make the medium 
more opaque. Instead the opacity saturates at a well de- 
fined value, when dipole-dipole interactions within the 
medium dominate the light-matter interactions. In ad- 
dition, when the layer thickness is less than A/4, we ob- 
serve a dramatic change in behaviour whereby the dipo- 
lar interaction doubles in strength, resulting in a differ- 
ent opacity that saturates at a value that is half that of 
the conventional case. This effect is independent of the 
transition dipole moment, and so is predicted to occur in 
other systems of two-level interacting dipoles. This result 
could have important consequences for the generation of 
non-classical light in confined geometries and in future 
work we will investigate the photon statistics associated 
with the transition between a 2D and 3D dipolar layer. 
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